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Abstract 



i-rt '. We study the influence of boundary conditions on the finite-size corrections of a one-dimensional 

Q I (ID) quantum spin model by exact and perturbative theoretic calculations. We obtain two new 



infinite sets of universal amplitude ratios for the finite-size correction terms of the ID quantum 



^j>. , spin model of A^ sites with free and antiperiodic boundary conditions. The results for the lowest 

o : 

two orders are in perfect agreement with a perturbative conformal field theory scenario proposed 



by Cardy [Nucl. Phys. B 270, 186 (1986)]. 
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I. INTRODUCTION 



Universality [3|-l9|, finite-size scaling |4|-|l6|. and finite-size corrections |17|-|29| for critical 
lattice systems have attracted much attention in recent decades. It has been found that 
critical systems can be classified into different universality classes so that the systems in 
the same class have the same set of critical exponents, whose values depends only on the 
global properties of the system such as spacial dimensions, number of components of the 



a. 



The hypothesis of 



order parameter, the range of interaction, and the symmetry group 

universality has much stronger implications and it is possible to show that models belonging 

to the same universality class also share the same set of universal finite-size scaling functions 



(UFSSFs) and amplitude ratios 



fl-llS 



21 



30 



3l|, whose values are independent of the 



microscopic structure of interactions. However, using a histogram Monte carlo simulation 
method [3^ , Hu, et al. confirmed for percolation the well-established dependence of finite- 



size sea 
|lO 



11 



ing functions on boundary conditions and on the shape of the domains of the systems 



14 



16| . Based on exact partition functions of the Ising model on square, triangular. 



and honeycomb lattices 33|, Wu, Hu and Izmailian have found similar results for the Ising 



34|. 



model 

In 1967-1969 Ferdinand and Fisher [17| calculated exact finite-size corrections (FSCs) 
for the free energy and its derivatives of the dimer and Ising models. In 1997, Ziff, Finch 



and Adamchik 



18| used Monte Carlo simulations to calculate FSCs for cluster numbers of 



two-dimensional random percolation models. In 2002, Caselle, et al. 19| used conformal 
field theory to study correction terms for the free energy and its derivatives of the Ising 
model. Based on the connections between lattice phase transition models and correlated 
percolation models [35|, in 1999 Hu, et al. |20| calculated FSCs for cluster numbers of the 
g-state Potts model for q being 1, 2, 3 and 4. In such studies, they extended the calculations 
of Ferdinand and Fisher for the Ising model on torus [17J to higher orders. Such studies 
inspire further research on FSCs for the Ising [21!] and the dimer models. 

Based on Kronecker's double series and exact asymptotic expansions for exact partition 



unctions 



26 



22| . exact finite-size corrections for the Ising model 22N25| and the dimer model 



29| on planar lattices with various boundary conditions have been obtained. Jt has 



been found that such correction terms also depend on the boundary conditions. In ^], we 
have found exact universal amplitude ratios for finite-size corrections of the two-dimensional 



Ising model on square, plane triangular and honeycomb lattices and a quantum spin chain on 
lattices with periodic boundary conditions. It is of interest to know how universal amplitude 
ratios depend on boundary conditions. In the present paper, we will address this question 
for a one dimensional quantum spin model with a varying parameter 7. We have calculated 
exact amplitude ratios for the quantum spin model with different boundary conditions and 
found that such ratios depend on boundary conditions and are independent of 7. 

The quantitative description of the universality classes of critical behavior is the main 
goal of quantum field theory when it is applied to statistical mechanics. In principle, the uni- 
versal quantities can be computed from the quantum field theory encoding the fundamental 
symmetries of the system. 

The universality class can describe the critical behavior of many physical systems, which 
undergo a second order phase transition. The criteria by which the critical systems can be 
classified into different universality classes is a problem of much academic interest. Two- 
dimensional (2D) critical systems are parameterized by the conformal anomaly c which 
is the central charge in the Virasoro algebra 36l-l38| . The conformal anomaly c can be 



obtained directly from the finite size corrections to the free energy for a 2D classical system 
on infinitely long but finitely wide strip at a conformal invariant critical point. 

The asymptotic finite-size scaling behavior of the critical free energy (Jn) and the inverse 
correlation lengths (^j~^) of an infinitely long 2D strip of finite width A^ at criticality has the 
form y, |4o| 

lim N\fM - /oo) - 2NUrf = A, (1) 

iV— !>oo 

hm iVC' = D,, (2) 

where foo is the bulk free energy, f^urf is the surface free energy and A and Di are the 
universal constants, but may depend on the boundary conditions (BCs). In some 2D ge- 



ometries, the values of A and Dj are known 



39N41| , to be related to the conformal anomaly 



number (c), the conformal weight of the ground state (A), and the scaling dimension of the 
2-th scaling field (xj) of the theory 

A = An(^-A), A = 2nx„ (3) 



.24 
for periodic or antiperiodic BCs and 



A = 7i{^-A), Di = 7ix„ (4) 



for free BCs. For the 2D Ising model, we have c = 1/2. The conformal weight of the ground 
state A, and the scahng dimension Xi are the universal quantities which depends on the 
BCs: for periodic BCs (A = 0,a;i = 1/8, X2 = 1), for antiperiodic BCs (A = 1/16, Xi = 
3/8,0:2 = 1), and for free BCs (A = O^i = 1/2, xs = 2). 



Quite recently, Izmailian and Hu J2l|] studied the finite size correction terms for the free 
energy per spin and the inverse correlation lengths of critical 2D Ising models on N ^ 00 
lattice and ID quantum Ising chain with periodic BCs. They obtained analytic expressions 
for the finite-size correction coefficients a^, hk and Ck in the expansions 

iV[/(iV)-/oo] = E^' (5) 



fc=i 






jY2fc-i • 
fc=i 

fc=i 
and find that although the finite-size correction coefficients a^, hk and Ck are not universal, 
the amplitude ratios for the coefficients of these series are universal and given by 

r^i^k) = h= f'~^ , (8) 

^ ' ttk 22'=-! -1 ^ ' 

'•'<*' = 5 = (2--'-l)B..' (^' 

where i?„ is the n-th Bernoulli number {B2 = 1/6, B4 = —1/30, . . . ). 

In the next section of the present paper we will present exact calculations for a set of 



universal amplitude ratios for the ID quantum spin model [42], which is the quantum version 
of the classical 2D Ising model, with free and antiperiodic BCs. We obtain analytic equations 
for Ofc, bk and Ck in the expansions given by Eqs. ([5]), IQ and (jTj) and find that universal 
amplitude ratios for the ID quantum Ising model with antiperiodic boundary condition are 
given by 

rsik) ^ ^ ^ ^^"-J\^--^^ (10) 

^ ^ ak 22fc-ifi2fc ^ ^ 

re{k) = ^ = -|^, (11) 

O-k -D2fc 



while for free BCs we obtain 




rs{k) = 


ak 


re{k) -- 


Ck 

ak 



1^ (12) 

4fc(3^^-^ + 1) . . 

As far as we know, no previous RG arguments, analytic calculations, or numerical studies 
predict the existence of this whole set of universal amplitude ratios. 

II. ONE-DIMENSIONAL QUANTUM SPIN CHAIN 

The ID quantum spin chain, the hamiltonian limit of the classical 2D Ising model (see 



43[ and references therein), belongs to the most frequently investigated and best understood 
systems in statistical physics. We will calculate universal amplitude ratios for this systems 
with anti-periodic boundary conditions and free boundary conditions. 

A. Anti-periodic boundary conditions 

Let us first consider the antiperiodic BCs. On a chain with A^ sites the Hamiltonian is 

given by 

\ ^ 1 ^ 

^ = -^ E < - II E [(1 + ^)<+i< + (1 - ^>i^i<] ' (14) 

^ n=l ^ p^ 

which was exactly solved by Katsura [42|. Here a^,a^ and a^ are the Pauli spin matrices. 



The antiperiodic BCs imposed on Eq. flT^ are 

cr%+i = -ol and a^_^^ = -o\. 



The phase diagram is well known [4J]. For all 7 (0 < 7 < 1), there is a critical point 
at Ac = 1, which falls into the 2D Ising universality class. For 7 = 1 it is also called ID 
transverse Ising model. Thus, by the introduction of a parameter 7 we could study different 
models in the same universality class. The Hamiltonian of Eq. flT^ can be diagonalized by 
a Jordan- Wigner transformation as 

E = Y,k^{rilr]^-\l'l), (15) 



where r]l,r]k are fermionic creation and annihilation operators and 



Afe = ■\/(cosA;-A)V7^ + sin2fc, (16) 

is the lattice dispersion relation. At the critical point Ac = 1 one then obtains A^ = 2ip {^) ■ 
Here 

tp{x) = Wsin^(x) ^ — sin^(x). (17) 

The critical ground-state energy, Eq , corresponds to the antiperiodic BC has the value 45 1 

E, /'Tcm\ v^ / o Trm 



E, 



(A) 



Tim 7^ — 1 . 4 7im 



m=0 m=0 



sm 



T 



N 



(18) 



The energy gaps As and Ae are given by 

7V-1 



....,,(_L),^,(^)_, 



m=0 '- 



/2m + 1 



2N 



--K 



A(^) = 2^|J 



n 



2 \ / sin 



2 TT 72 - 1 77 

7; — sm — . 

2 N 



(19) 



(20) 



.NJ \i N r 

Note, that the ground state energy Ei^, the first energy gap {E\ — Eq = A^) and the 
second energy gap (£'2 — Eq = Ag) of a quantum spin chain are, respectively, the quantum 
analogies of the free energy f{N), inverse spin-spin correlation length ^~^(A^), and inverse 
energy-energy correlation length ^~^(A^) for the Ising model; that is, 

Nf{N) ^ -Eo, C\N) ^ A„ and C\N) ^ A,. 

To write Eq , Ai , and Ai in the form of Eqs. (|5]), ([6]), and ([7]), we must evaluate 
Eqs. ( ITS]) . (IT^ . and ( 12U]) asymptotically. These sums can be handled by using the Euler- 
Maclaurin summation formula J46|. Suppose that F{x) together with its derivatives is 
continuous within the interval (a, b). Then the general Euler-Maclaurin summation formula 
states 

N-l , ^h -^ T/fc-l 



J2 F{a + nh + ah) = - F(r) dr + ^ -rrBfc(a) {F^'-'\b) - F^''-'\a)) 

n=0 "^" fc=l 



(21) 



where 0<a<l,/i = (6 — a)/A^ and Bfc(a) are so-called Bernoulli polynomials defined in 
terms of the Bernoulli numbers Bp by 

k 



B,(a) = 5^B, 



p=0 



fc! 
(/c — p)\p\ 



a 



k—p 



(22) 



Indeed, -Bn(O) = i?„- Bernoulli polynomials satisfy the identity: 

5„(l/2) = (2i-"-l)5„ (23) 



By expanding the exact solution of Eq. ( 1T^ . Henkel J45| has obtained several finite-size 
correction terms to the ground-state energy Eq and to the first energy gap Ei — Eq. We 
have extended the calculations to arbitrary order and found that 



-r>.A^"o^E^(^) ■^- 



7r_l/l_4\/7r\3 

m'TEW's) \2n) + ■ ■ ■ ' ^^^ 

2fe _ n / ^ X 2k-i 



A(^) = V 4fc - 2i?,.(2^^ - 1) / ^ y- -(,,_,) 
^ (2fc)! V2Ary ^ 



fc=i 

Svr 9/1 4\ / TT \3 

4iv + 8(y-3J(2iv) +■■■• (2«' 

4A; /7r\2fe-i 
:^ (2A0! liV. 

27r /I 4\ /7r\3 
where ■0(2'^"^) = {d?^~^il){x) / dx^'^~^^ and ao is an non-universal number 






ao = — / tp{x)dx = 1 + arccos7/(7V 1 ^ 7^) /^- (27) 

The ratios of the amplitudes of the A^^^^'^^^) correction terms in the spin-spin correlation 
length, energy-energy correlation length and the free energy expansion, i.e. bk/dk and Ck/ak, 
are 7- independent and are given by Eqs. flTOj) and (ITTl) . respectively, thus confirming the 
universality of this ratios. 

B. Free boundary conditions 

Let us now consider the transverse Ising model on a ID lattice of A^ sites with free BC 
and with the Hamiltonian 

Af N-1 

H = -XY,<-2Y.«+v (28) 

n=l n=l 



In the limit N ^ oo, the ground state is singular at A = 1. The Ising chain in a transverse 
field has been studied in great detail (see e.g. |42|, |47|). The Hamiltonian of Eq. (I28l) may 
be re-expressed in the diagonal form 

H = J2^k{vivk-l/2) (29) 

k 

with dispersion relation 



Ak = ^{X-l)^ + 4Xsm\k/2). (30) 

or free BCs one finds that the allowed values of k are determined by the secular equation 



47| 



A-i = sin[(A^ + l)k]/ sm{Nk). (31) 

At the critical field Ac = 1, the secular equation reduced to 

tan(fcA^) = cot(fc/2), (32) 

so that 

(2m+2h ^ = o,l,...,iV-l. 

2N + 1 ' ' ' ' 

and for the dispersion relation one then obtain A^ = 2^? f m+i ) ' where 

(p{x) = sin(x). (33) 

The ground-state energy Eq , the energy gaps Ai and Ai are given by 47 1 



E^ = -y^ ( ^^^^TT ) = i ( 1 - cosec^^ ) , (34) 



Af) = 2^(^^ =2 sin 



n 



2{2N+1)' 



(35) 



^'^''-'i--2i2^)^''''m^)^- ^''^ 



The asymptotic expansion of these quantities can be written in the following form 

)2fe-l / \ 2fc-l 



k=l 



B2k{l 



IT 



{2k)\ 

7 



TT 



Af ' = E 



24(2iV+l) 
2 



n 



^ {2k- 1)! V4A^ + 2 



90 V8A^ + 4 

2fc-l 



AN + 2, 

3 

+ ... 



V 



(2fc-l) 



(37) 



TT 



1 



2N+1 3 V 4iV + 2 



vr 



Af ' = E 



2(3 



2fc-l 



1) 



k=l 

in 



{2k-l)\ 

28 
3" 



TT 



+ ..., 
2fc-l 



(2k-l) 



(38) 



4Ar + 2 

3 

TT 



^ 



(2fc-l) 



4A^ 



+ 



(39) 



2A^ + 1 

where, yj^^'^"^) = ((i^'^~^v^(x)/(ix^'''"^)^^g, (p{x) = sinx, and ao = ^ J^ (p{x)dx = 2/tt. 

Equations (137|) . ( 138|) . and ( 139|) imply that the ratios of the amplitudes of the A^^^^'^^i) 
correction terms in the spin-spin correlation length, the energy-energy correlation lengths, 
and the free energy expansion, i.e. bk/a^ and Ck/ak, should not depend in detail on the 
dispersion relation as given by Eqs. ( TT2l) and ( IT3l) . 

The leading terms of Eqs. (EH) - ([261), (EZD - ([39!) are consistent with Eqs. (^i) - (^, i.e. 
Oi, 6i and Ci are universal. Equations ([1]) and ([2]) implies immediately that their ratio is 
also universal, namely rs{l) = Di/A and re(l) = D2/A, which is consistent with Eqs. ([8]), 
(|T0|), and Eq. ([12]) for the case k = l 



rs{l) 



£1 
A 



and with Eqs. ([9]), ([TT]) . and Eq. ([T3|) for the case A; = 1 



3 


for periodic BC 


-9/2 


for antiperiodic BC 


24 


for free BC 



(40) 



re(l) = ^ 



24 



-12 



96 



for periodic BC 
for antiperiodic BC 
for free BC 



(41) 



III. PERTURBATIVE CONFORMAL FIELD THEORY 



.N/2 

H = Hc + 2_,9p I 4>p{v)dv, 

J-N/2 



(42) 



The finite-size corrections to Eos. fl) and ([2]) can be calculated by the means of a 
perturbative conformal field theory 48|, |49|]. In general, any lattice Hamiltonian will contain 
correction terms to the critical Hamiltonian H^ 

fAf/2 
p ~'-N/2 

where Qp is a non- universal constant and 0p(f) is a perturbative conformal field. Below 
we will consider the case with only one perturbative conformal field, say (piiv). Then the 
eigenvalues of H are 

-N/2 

(43) 



En = En,c + 9il < n\(l)i{v)\n > dv + 

'-N/2 



where En.c are the critical eigenvalues of H. The matrix element < n\(j)i{v)\n > can be com- 
muted in terms of the universal structure constants (Cnin) of the operator product expansion 
48| : < n\(j)i{v)\n >= {2tt/N)^' Cnini where Xi is the scaling dimension of the conformal field 
(t>i{v). The energy gaps (A„ = En — Eq) and the ground-state energy {Eq) can be written as 



—Xn + 2'Kgi[Cnln - Cm) I -JT 



Xl-l 



+ ... 



Eq 



2n 
Eo,c + ^ngiCoio ( — 



Xl-l 



+ .... 



(44) 
(45) 



For the 2D Ising model, one finds 21[ that at least two (and probably infinitely many) 
perturbative conformal fields are necessary to generate all finite-size corrections terms. Nev- 
ertheless, the leading finite-size corrections (1/A^^) can be described by the Hamiltonian 
given by Eq. (H2|) with a single perturbative conformal field (f)i{v) 



LlJv) 



'-2 



L^gl"^) '^ith 



scaling dimension x/ = 4 50 1. 

In order to obtain the corrections we need the matrix elements < n\E^_2{v) + L?_2(f )|n >, 
which have already been computed by Reinicke 51| : 



< A + r|L^2l^ + r > 



< r\L_r\r > 



27r 

271 



49 +(^ + ^)^ 5 , r(2A + r)(5A + l) 



11520 
49 



24 (A+1)(2A + 1^ 



49 
+ Z^r(2r2-3) 



;46) 

(47) 



11520 120 
The universal structure constants C212, Ciii and Cqiq can be obtained from the matrix element 

< n\L^_,{v) + L^_,{v)\n >= {27^/Nf' C^m. (48) 



10 



where xi = 4 is the scahng dimension of the conformal field L'^_2{v) + L'^_2{v). 

At the critical point Ac = 1 the spectra of the Hamiltonian f lT4|) with periodic and 
antiperiodic BC are built by the irreducible re pre sentation A, A of two commuting Virasoro 
algebras L„ and L„ with central charge c = ^ [41]. We denote by A the highest weight, and 
by A + r, the r-th level having degeneracy (i(A, r) of one irreducible representation of the 
Virasoro algebra. A state will be labelled by \n >~ |A + r, A + f>. The possible values of 
A, A are (0, 0), (|, |), (^, ^) for periodic and (0, |), (|, 0), (^, ^) for antiperiodic BCs. In 
the case of free BCs the spectra can be understood in terms of irreducible representations 
A of a single Virasoro algebra with possible values of A are 0, |. 

The ground state |0 >, first excited stat e | 1 >, and second excited state |2 > depends on 



the boundary conditions and given by 48 



5J: 



for periodic BCs; 



|0 > = I A = 0, r = 0; A = 0, f = >, (49) 

|1> = |A = ^,r = 0;A = ^,f = 0>, (50) 

|2> = |A = i,r = 0;A = i,r = 0>, (51) 



|0> = |A = -^,r = 0;A = -^,f = 0>, (52) 
lo lo 

|1 > = |A = 0,r = 0;A = ^,f = 0>, (53) 

1 - 1 

|2 > = I A = — , r = 0; A = — , f = 1 >, (54) 

' ' 16 16 ^ ' 



for antiperiodic BCs; and 



|0> = |A = 0,r = 0>, (55) 

|1> = |A = l,r = 0>, (56) 

|2> = |A = 0,r = 2> (57) 



for free BCs. 

After reaching this point, one can easily compute the universal structure constants Coio, 
Cm and C212 for all three boundary conditions. The values of Cqzo, Cm, C212 can be obtained 
from Eqs. P6|) - fl571) and given by: 



Cozo = 49/5760, Cm = -7/720, C2Z2 = 1729/5760, (58) 

11 



for periodic BCs; 



a 



010 



-7/720, Cm = 889/5760, C212 = 833/720, 



for antiperiodic BCs; and 



Cqio = 49/11520, Cm = 1729/11520, C212 = 47089/11520, 



(59) 



(60) 



for free BCs. Equations (jHj) and (H5l) implies that the ratios of first-order corrections 
amphtudes for (A„) and {—Eq) is universal and equal to {Coio—Cnin)/CoiQ, which is consistent 
with Eqs. ([H]), dlQ]), and Eq. (^2l for the case k = 2 



rsi2) 



Cow ~ Cm 



c 



010 



15/7 for periodic BCs, 

135/8 for antiperiodic BCs, 
-240/7 for free BCs 



and with Eqs. ([9]), ([II]), and Eq. (^ for the case k = 2 



re{2) 



Cqio — C212 



C, 



010 



-240/7 for periodic BCs, 
120 for antiperiodic BCs, 

-960 for free BCs. 



(61) 



(62) 



IV. CONCLUSION 



In this paper we discuss the influence of the boundary conditions on the finite-size correc- 
tions of the 2D Ising model in the extreme anisotropic or quantum-Hamiltonian limit. We 
have calculated various universal amplitude ratios and find that such result are in perfect 
agreement with a perturbated conformal field theory scenario proposed by Cardy 48 1. 

The results of this paper inspire several problems for further studies: (i) Further work 
has to be done to possibly evaluate exactly all finite-size correction terms from perturbative 
conformal field theory, (ii) Can one obtain from the perturbated conformal field theory 
the value of the universal amplitude ratios {rs{k) and re{k)) for k > 27 (iii) How do such 
amplitudes behave in other models, for example in the three-state Potts model? (iv) Our 
results also present new challenges to scientists working on numerical studies of critical 
phenomena. For example, it is of interest to present accurate numerical evidences about 
whether the Ising model on a two-dimensional lattice with crossing bonds has the same set 
of amplitude ratios. 



12 
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